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ABSTRACT 

Context. Galaxy merging is widely accepted to be a key driving factor in galaxy formation and evolution, while the feedback from 
actively accreting nuclei is thought to regulate the black hole-bulge coevolution and the star formation process. 

Aims. In this context, we focused on ISXPS J050819.8+172149, a local (z = 0.0175) Seyfert 1.9 galaxy (Lbol ~ 4 x 10"^^ ergs s“*). 
The source belongs to an infrared-luminous interacting pair of galaxies, characterized by a luminosity for the whole system (due to the 
combination of star formation and accretion) of log(LiR/Lo) = 11.2. We present here the first detailed description of the 0.3 - lOkeV 
spectrum of ISXPS J050819.8+172149, monitored by Swift with 9 pointings performed in less than 1 month. 

Methods. The X-ray emission of ISXPS 1050819.8+172149 is analysed by combining all the Swift pointings, for a total of ~ 72ks 
XRT net exposure. The averaged Swift-BAT spectrum from the 70-month survey is also analysed. 

Results. The slope of the continuum is F ~ 1.8, with an intrinsic column density of ~ 2.4 x 10^^ cm“^, and a de-absorbed luminosity 
of ~ 4 X 10**^ ergs s“* in the 2 - 10 keV band. Our observations provide a tentative (2. Icr) detection of a blue-shifted Fe xxvi absorption 
line (rest-frame E ~ 7.8keV), thus suggesting the discovery for a new candidate powerful wind in ISXPS J050819.8+172149. The 
physical properties of the outflow cannot be firmly assessed, due to the low statistics of the spectrum and to the observed energy of the 
line, too close to the higher boundary of the Swift-XRT bandpass. However, our analysis suggests that, if the detection is confirmed, 
the line could be associated with a high-velocity (uqui ^ 0.1c) outflow most likely launched within 80 rj. To our knowledge this is 
the first detection of a previously unknown ultrafast wind with Swift. The high column density suggested by the observed equivalent 

width of the line (EW-230 eV, although with large uncertainties), would imply a kinetic output strong enough to be comparable to 

the AGN bolometric luminosity. 
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1. Introduction 


The observational evidence for the presence of inactive Super 
Massive Black Holes (SMBHs; Mbh ~ 10® - 10*° Mq) at the 
centre of most, if not all, the local galaxies, and the observed 
correlation between several properties of the galaxy’s bulge and 
the central SMBH mass (Ferrarese & Merritt 2000; Gebhardt 
St al. 2000), suggest that the SMBH accretion and the assembly 
of the galaxies bulges are intimately related (see Kormendy & 
Ho 2013, for a recent review). Funnelling of gas in the nuclear 
regions, as triggered by galaxy interactions, can activate both 
efficient accretion onto the SMBH, and a burst of star forma¬ 
tion. A key ingredient in regulating their evolution should be the 
feedback from the Active Galactic Nuclei (AGN); being conser¬ 
vative, while building its mass the SMBH can release an amount 
of energy larger than ~ 30 times the binding energy of the host 
bulge (see a review in Fabian 2012). Even if only a small frac¬ 
tion of this energy is transferred to the gas in the galaxy, then 
an active nucleus can have a profound effect on the evolution 
of its host (Di Matteo et al. 2005). Powerful (kinetical or radia¬ 
tive) outflows of gas driven by luminous quasars are invoked as 


* Based on observations obtained with the Swift satellite. 
** E-mail: lucia.ballo@brera.inaf.it (LB) 


a key mechanism to blow away the gas in the galaxy and thereby 
quench star formation, coincidentally starving the SMBH of fuel 
(King & Pounds 2015). 

AGN winds with a range of physical properties have been 
revealed by observations at various energies, from radio up to 
X-rays. Outflows of molecular or neutral atomic gas, with veloc¬ 
ities up to ~ 1000 - 2000 km s * and extending on kpc scales, 
have been observed at mm (e.g., Feruglio et al. 2010; Cicone 
et al. 2014) and radio frequencies (e.g., Morganti et al. 2005; 
Teng et al. 2013) in a few dozen AGN in dusty star form¬ 
ing sources and/or radio galaxies. Mass outflows of ionized gas 
with similar velocities at distances consistent with the narrow 
line region zone have been detected in the optical/ultraviolet 
(UV), both in the [Oiii] emission line profiles (e.g., Crenshaw 
& Kraemer 2005; Cano-Dlaz et al. 2012; Cresci et al. 2015), 
and through the observation of broad absorption line systems 
(e.g., Dai et al. 2008; Borguet et al. 2013). In X-rays, mildly 
ionized warm absorbers are observed in more than half of unob¬ 
scured AGN (Crenshaw et al. 2003; Blustin et al. 2005). The ob¬ 
served velocities of ~ 500 - 1000 km s * imply a kinetic power 
rather low when compared to the bolometric luminosity. How¬ 
ever, Crenshaw & Kraemer (2012) found that, when summed 
over all the absorbers, the total power carried by these structures 
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can reach, in some cases, the minimum level required for AGN 
feedback (L/Lboi ~ 0.5 - 5%; e.g., Di Matteo et al. 2005; Hop¬ 
kins & Elvis 2010). In the last years, highly blue-shifted Fe K- 
shell absorption lines at rest-frame energies E >1 keV, observed 
in XMM-Aewton or Suzaku spectra of luminous AGN, revealed 
the presence of high column density (Ah ~ 10^^ cm“^) and fast 
(v > 0.1c) winds (e.g., Tombesi et al. 2010; Gofford et al. 2013; 
Tombesi et al. 2015). Their derived kinetic power is systemati¬ 
cally higher than the minimum fraction of bolometric luminos¬ 
ity required by AGN feedback models in order to regulate the 
growth of an SMBH and its galactic bulge (e.g., Hopkins & Elvis 
2010). Only very recently, outflows over a range of scales have 
been detected and studied within the same source, thus allowing 
us to explore the connection between large-scale molecular out¬ 
flows and accretion-disk activity (Feruglio et al. 2015; Tombesi 
et al. 2015). 

Finding observational evidence of the effects of such 
accretion-related feedback and characterising the magnitude of 
mass outflows from AGN are among the major challenges of 
the current extragalactic astronomy. Here we present 9 coadded 
Swift observations (performed in less than 1 month) of the in¬ 
teracting infrared (IR) galaxy ISXPS J050819.8H-172149. The 
spectrum, the first for this source covering the energy range 
E ~ 1 - lOkeV, provides us with a tentative detection of a pos¬ 
sible new ultrafast wind (uout ~ 0.1c). To our knowledge, this is 
the first time that a previously unknown outflow is revealed by 
Swift. 

The source is described in Sect. 2, while the analysis of 
the Swift-XRT data is presented in Sect. 3 and the results are 
summarised in Sect. 4. Throughout the paper we assume a flat 
ACDM cosmology with Hq - 71 km s“* Mpc“\ Qa = 0.7 and 
Gm = 0.3. 

2. 1SXPS J050819.8+172149 

The position of the Swift source ISXPS J050819.8+172149 
(Evans et al. 2014) is coincident with the center of CGCG 468- 
002 NEDOl, a local (z - 0.0175) galaxy belonging to the in¬ 
teracting pair CGCG 468-002 (projected distance between the 
center of the galaxies of ~ 29.4", corresponding to ~ lOkpc at 
the source redshift). 

The high 8 - lOOOyum luminosity of the system, 
log(LiR/Lo) = 11-2 (Armus et al. 2009), implies a classifi¬ 
cation as Luminous Infrared Galaxy (LIRG, defined as hav¬ 
ing Lir ^ lO" '^Lq; Sanders & Mirabel 1996). Included for 
this reason in the Great Observatories All-Sky LIRG Survey 
(GOALS; Armus et al. 2009), the system has been targeted with 
several observational facilities, and multiwavelength informa¬ 
tion (both spectroscopic and photometric) has been collected. 
The available data span from the UV (GALEX; Howell et al. 
2010) up to the mid-IR (Spitzer/IRS, IRAC, and MIPS; Diaz- 
Santos et al. 2010, 2011; Petrie et al. 2011; Inami et al. 2013; 
Stierwalt et al. 2013, 2014, see also Valiante et al. 2009; Alonso- 
Herrero et al. 2012,2013) and far-IR (HERSCHEL/PACS; Diaz- 
Santos et al. 2013) bands. Photometry from the 2MASS Redshift 
Survey is reported by Huchra et al. (2012). At radio wavelengths, 
the relatively strong (S y ~ 34mJy) NVSS emission (Condon 
et al. 1998), detected halfway between the nuclei and slightly 
elongated in their direction, is probably due to the combined con¬ 
tribution of both galaxies. 

The system is an early-stage merger (Stierwalt et al. 2013) 
known to host an AGN optically classified as Seyfert 1.9 (see 
e.g. Motch et al. 1998; Veron-Cetty & Veron 2001; Kollatschny 


Table 1. Swifl-XKT monitoring observation log for 
ISXPS J050819.8+172149. Net count rates are in the 0.3 - lOkeV 
energy range. 


Obs. ID 

Start date 

Net count rate 
[10^^ counts s“*] 

Net exp. time 

[s] 

00049706003 

2014-10-15 

7.5 + 0.4 

4235 

00049706004 

2014-10-16 

13.1+0.3 

14740 

00049706005 

2014-10-21 

10.0 + 0.7 

2308 

00049706006 

2014-10-22 

7.3 +0.3 

7325 

00049706007 

2014-10-23 

7.5 +0.3 

6533 

00049706008 

2014-10-27 

5.8 +0.3 

6478 

00049706009 

2014-10-28 

5.9+ 0.3 

6218 

00049706010 

2014-10-30 

4.6+ 0.3 

6333 

00049706011 

2014-11-11 

7.6+ 0.2 

17950 


et al. 2008). Mainly based on the detection of the high-ionization 
line [Ne v] in the SpitzerflRS spectra, the Seyfert 1.9 nucleus has 
been associated with the western galaxy of the system, coinci¬ 
dent with ISXPS J050819.8-H172149 (Petrie et al. 2011; Alonso- 
Herrero et al. 2012; Stierwalt et al. 2013, 2014, see also the op¬ 
tical classification reported by Alonso-Herrero et al. 2013). The 
Spitzer spectra are suggestive of a relatively unobscured AGN 
(consistent with the optical classification as Seyfert 1.9), ener¬ 
getically important in the mid-IR (Stierwalt et al. 2013). From 
a decomposition of the mid-IR spectra into AGN and starburst 
components, Alonso-Herrero et al. (2012) estimated a bolomet¬ 
ric luminosity due to the accretion of Lboi ~ lO^^L©. Instead, 
no signature of active accretion is found for the eastern galaxy, 
which shows all the typical properties of a star forming source 
both in the mid-IR spectra and from the UV photometry. 

The BH mass in ISXPS J050819.8+172149 is among the 
highest observed in local LIRGs, Mbh ~ 1-15 x 10^ Mq 
( as calculated from the velocity dispersion of the core of the 
[Oiii]/15007 line; Alonso-Herrero et al. 2013). This implies that 
the black hole is radiating at a low fraction of its Eddington lumi¬ 
nosity', with Eddington ratio logTedd = logTboi/TEdd-2.5. 

The star formation rate (SFR) derived for this source falls in 
the lower tail of the distribution found for local LIRGs (nuclear 
SFR ~ 1 - 3 Mq yr“', integrated SFR ~ 3 - 4 Mq yr“*, Alonso- 
Herrero et al. 2013), implying a ratio between SFR and BH ac¬ 
cretion rate of log SFR jmsH ~ 2 (by assuming a mass-energy 
conversion efficiency e = 0.1), similar to the values found for 
Seyfert galaxies. 

3. XRT data analysis 

Besides the basic analysis reported in the Swift-XRT point 
source catalogue (Evans et al. 2014), in the soft-medium X-ray 
energy range (E < lOkeV) the only published information up 
to now comes from the ROSAT All Sky Survey (see Kollatschny 
et al. 2008, and references therein), providing a soft X-ray lu¬ 
minosity (de-absorbed by our Galaxy) of To.i-2.4keV ~ 6.6 x 
10'"ergss-'. At higher energies, ISXPS J050819.8+172149 
has been detected by the Swift-Burst Alert Telescope (BAT; 
see Baumgartner et al. 2013), while Ackermann et al. (2012) 
reported only a Fermi 95% confidence level upper limit of 
To.i-iOGeV ~ 1-6 X lO^'^ergs s^'. 

' The Eddington luminosity, defined as Ledd = 47rGMBHfnpC/cr-j- 
1.3 X Iff* MbuIMq [ergs s“'], represents the exact balance between in¬ 
ward gravitational force and outward radiation force acting on the gas, 
assumed to be of ionised hydrogen in a spherical configuration. 
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Recently, our group has been awarded a Swift (Gehrels et al. 
2004) program for this source (PI P. Severgnini): 9 pointings per¬ 
formed with the X-ray telescope (XRT; Burrows et al. 2005) in 
the standard photon counting (PC) mode between 2014-10-15 
and 2014-11-11, for a total of ~ 72 ks net exposure (ObsID from 
00049706003 to 00049706011; see Table 1). 

We generated images, light curves, and spectra, including the 
background and ancillary response files, with the online XRT 
data product generator^ (Evans et al. 2007, 2009); the appropri¬ 
ate spectral response files have been identified in the calibration 
database. The source appears point-like in the XRT image and 
centered at the position of the western nucleus, without any ev¬ 
ident elongation toward the position of the second galaxy. We 
note that at the angular resolution of XRT, 18" half-power di¬ 
ameter, the emission of the two galaxies, located at a distance of 
~ 29.4", can be resolved; by assuming a power-law model with 
r = 2, we estimated a 3cr upper limit to the 0.3-10keVemission 
at the eastern source position of ~ 7 x 10“'"^ ergs cm“^ s"'. 

Source events were extracted from a circular region with a 
radius of 20 pixels (which corresponds to an Encircled Energy 
Eraction of 90%, Moretti et al. 2005, 1 pixel ~ 2.36 arcsec), 
while background events were extracted from an annulus re¬ 
gion centred on the source with inner and outer radii of 60 and 
180 pixels, respectively; all sources identified in the image were 
removed from the background region. 

We extracted a light curve binned at the duration of each 
individual observation. ISXPS J050819.8+172149 was detected 
in all observations, with 0.3-10keV signal-to-noise ratios {SIN) 
ranging from 15 to 44. The average count rates in the total (0.3 - 
lOkeV), soft (0.3 - 2keV), and hard (2 - lOkeV) XRT energy 
ranges are ~ 0.09, ~ 0.03, and ~ 0.07 counts s“\ respectively. 
Small deviations from these values, of a factor lower than 2, are 
observed in the light curves; however, there is no evidence of 
spectral variability in the ratio of count rates observed in the hard 
and soft bands. We do not find any significant pile-up problem. 

In order to increase the statistics, we co-added the XRT 
datasets. Source and background spectra were extracted from 
the merged event lists, and the former was binned in order 
to have at least 20 total counts per energy channel. The net 
count rates in the 0.3 - lOkeV, 0.3 - 2keV, and 2 - lOkeV 
energy ranges are (8.2 ± 0.1) x 10“^, (2.3 ± 0.1) x 10“^, and 
(5.8 + 0.1) X 10“^ counts s“\ respectively. The S/N achieved in 
the same energy ranges are 76, 40, and 64, respectively. 

Spectral fits were performed in the 0.3 -10 keV energy range 
using the X-ray spectral fitting package XSPEC (Arnaud 1996) 
vl2.8.2. Uncertainties are quoted at the 90% confidence level 
for one parameter of interest — 2.71). All the models dis¬ 
cussed in the following assume Galactic absorption with a col¬ 
umn density of A?H.Gai = 1.84 x 10^^ cm“^ (Kalberla et al. 2005). 
To model both Galactic and intrinsic absorptions we used the 
(z)pHABS model in XSPEC, adopting cross-sections and abun¬ 
dances of Wilms et al. (2000). 

A fit with a simple absorbed power law clearly provides a 
poor representation of the XRT data (;t'^/d.o.f. = 280.8/229); the 
photon index is T = 1.49 ± 0.09 and the column density is Ah = 
(1.8 + 0.2) X 10^^ cm“^. Above ~ 5 keV, i.e. in the energy range 
where the iron K complex is expected, residuals are present both 
in emission and in absorption (see Eig. 1), while a big bump 
below 1 keV suggests the emerging of an additional component 
in the soft band. Assuming that the neutral absorber only partial 
covers the central source (zpcfabs model in XSPEC) still leaves 
residuals in the soft band (x^lA.o.f. - 271.0/228); we found a 

^ http://www.swift.ac.uk/user_objects 


XRT: Abs. PL model 



Fig. 1. Residuals, plotted in terms of sigmas with error bars of size one, 
for the XRT spectrum, modelled with an absorbed power-law compo¬ 
nent (r ~ 1.5, Ah ~ 1.8 X lO^^cm^^). A deep absorption trough is 
seen in the iron K band, possibly associated with a broad absorption 
feature (dashed red line). The dotted red line marks the position where 
the neutral Fe Kcr emission line is expected (and marginally detected). 


covering fraction > 0.95, with T = 1.48 + 0.04 and Ah = (2.1 + 
0.1) X lO^^cm^^. 

In principle, the emission observed at low energies can be 
associated with the accreting nucleus (i.e., due to scattering off 
optically thin ionised gas) and/or the host galaxy (a soft ther¬ 
mal emission is a characteristic signature in all known starburst 
galaxies). Phenomenologically, the observed soft excess can be 
accounted for equally well by adding to the previous model: 
7) an unabsorbed power law with photon index tied to the pri¬ 
mary one (x^lA.o.f. - 258.4/228); the best-fit parameters are: 
r = 1.57 + 0.09, Ah = (2.2 + 0.2) x 10^2 cm^^, and ratio of 
the power-law normalizations A scatt/?l inn ~ 3%; or 2) a thermal 
component (mekal model in XSPEC; ;^f2/d.o.f. = 259.3/227): in 
this case we found r = 1.53+0.09, Ah = (1.9 ±0.2) x 10^2 cm“2, 
and kT - 0.3UQQgkeV. In the latter parametrisation, the lu¬ 
minosity attributed to the thermal component is T™j 2 kev ~ 
2.6xlO"^**ergs s“', which would imply a 57^7? ~ 0.5-2.9Moyr“' 
(e.g., Ranalli et al. 2003; Mas-Hesse et al. 2008; Mineo et al. 
2012), consistent with the nuclear star formation properties of 
the host (see Sect. 2). In both cases, the luminosity observed in 
the range covered by ROSAT, To.i-2.4keV ~ 5.6 x 10^' ergs s“', is 
in agreement with the value derived by Kollatschny et al. (2008). 

While the most plausible hypothesis is that both components 
contribute to the emission observed at low energies in the XRT 
spectrum, the quality of the present data does not allow us to dis¬ 
criminate between the two contributions, and even less to disen¬ 
tangle them. In the following, we assume an unabsorbed power 
law, checking that the inclusion of a thermal component in place 
of the power law does not affect the main results presented here. 

The addition of a narrow (50 eV) Gaussian emission line 
to the absorbed plus unabsorbed power laws results in an im¬ 
provement in the fit {Ax^ jAA.o.f. - 11.3/2). The line param¬ 
eters are E - 6.40 + 0.09 keV, consistent with neutral Ee Ka, 
7 = 1.5 ± 0.7 X 10“^ photons cm 2 s ', and EW~ 200eV. The 
strength of the Fe Ka line, coupled with the hard photon in¬ 
dex (T ~ 1.62), could suggest the presence of neutral reflec¬ 
tion (Reynolds et al. 1994; Matt et al. 1996, 2000). This pos¬ 
sible component was then included in the model by replacing 
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Table 2. Summary of the Swift (XRT+BAT) parameters of the best-fit model described in Sect. 3, where a high-energy absorption feature is 
superimposed to a continuum composed by an absorbed power law plus a reflection component, with the addition of a scattered power law. 


Direct, reflected and scattered continuum Absorption line 

1% T R % scan. E Ai EW logF2_iokeV log^2-i0keV A^d-Q.f. 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 


2.4 + 0.3 1.8 + 0.1 1.8!“-^ X 10“^ 1.4 + 0.4 2.1 +0.8 7.8 ± 0.3 < 270 -(1.3!^;^) x lO’^ -(230!]®“) -11.19 + 0.22 42.61 +0.18 241.5/228 


Notes. Errors are quoted at the 90% confidence level for 1 parameter of interest (Ax^ = 2.71). Col. (1): Intrinsic column density, in units of 
10^^ cm“^. Col. (2): Power-laws (absorbed and unabsorbed) and reflection component photon index. Col. (3): Absorbed power-law normalisation, 
in units of photons keV“* cm“^ s“*. Col. (4): Reflection fraction. Col. (5): Percentage of scattering fraction, defined as the ratio of the unabsorbed 
and absorbed power-law normalizations. Col. (6): Rest-frame energy centroid of the Gaussian absorption line, in units of keV. Col. (7): Absorption 
line width, in units of eV. Col. (8): Gaussian absorption line normalisation, in units of photons cm“^ s“'. Col. (9): Absorption line equivalent width, 
in units of eV. Col. (10): Observed flux (de-absorbed by our Galaxy) in the 2 - lOkeV energy band, in units of ergs cm“^ s“'. Col. (11): Intrinsic 
luminosity in the 2 - lOkeV energy band, in units of ergs s“*. Col. (12): and number of degree of freedom. 


the narrow Gaussian line with a pexmon component (Nandra 
et al. 2007), an additive model self-consistently incorporating 
the Compton-reflected continuum from a neutral slab combined 
with emission from Fe Ka, Fe Ky6, Ni Ka and the Fe Ka Comp¬ 
ton shoulder. 

To improve the determination of the slope of the pri¬ 
mary power law and the amount of reflection, it is funda¬ 
mental to know the shape and intensity of the emission at 
energies higher than lOkeV. Therefore we fitted the XRT 
data simultaneously with the averaged 5w//f-BAT spectrum of 
ISXPS J050819.8+172149 obtained from the 70-month survey 
archive^ (SWIFT J0508.1 + 1727). The data reduction and ex¬ 
traction procedure of the eight-channel spectrum is described in 
Baumgartner et al. (2013). To fit the pre-processed, background- 
subtracted BAT spectrum, we used the latest calibration re¬ 
sponse as of 2013 May. ISXPS J050819.8+172149 was detected 
in the 14 - 100keV band with a count rate of (3.9 ± 0.3) x 
10“^counts s ', which corresponds to a 14 - 195keV flux of 
2.6 ± 0.5 X 10“" ergs cm“^ s“' (Baumgartner et al. 2013). 

During the fit, the only free parameter of the pexmon compo¬ 
nent was the reflection scaling factor, R. We tied the pexmon pho¬ 
ton index and normalization to that of the primary power law, and 
we fixed the cutoff energy at 1000 keV (i.e., consistent with no 
measurable cut-off), the inclination angle at 60°, and the abun¬ 
dances of heavy elements at their Solar values. We also allow to 
vary the cross-normalisation factor between the XRT data and 
the average BAT spectrum. 

The baseline model then consists of an absorbed power law 
and an unabsorbed one plus a reflection component: the pho¬ 
ton index and intrinsic absorption of the primary emission are 
F = 1.75 + 0.09 (in agreement with the values typically observed 
in unobscured AGN; Piconcelli et al. 2005; Mateos et al. 2010; 
CoiTaletal. 2011)andAH = (2.4+0.2)x 10^^ cm“^, while the re¬ 
flection fraction and the strength of the scattered component are 
R ~ 1.1 and AscattMintr ~ 2%. The cross-normalisation factor 
between the XRT data and the average BAT spectrum is 0.96. 

Although this model provides an acceptable description 
of the broad-band (0.3 - 100 keV) X-ray continuum of 
ISXPS J050819.8+172149 ix^lA.o.f. = 250.2/231), it is not 
able to account for the residuals observed above 7 keV, confirm¬ 
ing the presence of a deep absorption trough at about 8 keV (see 
Fig. 2). Considering only the 41 XRT bins between 5 and 9keV, 
the fit statistics provided by the model is ;^f^/d.o.f. = 44.2/36. 
These residuals can be accounted for by adding a Gaussian ab¬ 
sorption line, with a rest-frame energy of £ = 7.8 + 0.3 keV and 


^ http://swift.gsfc.nasa.gOv/results/bs70mon/SWIFT_J0508.lpl727 


Baseline model: Abs. + unabs. PLs + reflected cont. 



j_I_I_I_I_I_I_I_I_ j _I_I_I_^_I_I_I_I_I_I_ 

6 7 8 9 10 

Rest Energy (keV) 


Fig. 2. Residuals, plotted in terms of sigmas with error bars of size one, 
for the XRT spectrum above 5.3 keV. The adopted model is an absorbed 
power law and a reflection component (F ~ 1.8, Ah ~ 2.4 x 10^^ cm“^, 
R ~ 1.1). An unabsorbed power law with a ~ 2% of intensity of the 
intrinsic one accounts for the emission below ~ 1 keV (actually, below 
the reported energy interval). The centroid of the absorption line when 
one more Gaussian is included in the model is marked with a red vertical 
line (absorption line at £ ~ 7.8 keV). 

an equivalent width EW= -(2301^^°) eV {Ax^lAA.o.f. = 8.3/3 
in the 5 - 9 keV interval). The line is marginally resolved; allow¬ 
ing its width to vary, we can set only an upper limit cte < 270 eV. 

The best-fitting parameters are reported in Table 2, while the 
unfolded XRT and BAT spectra are shown in Fig. 3. We note 
that the significance of the detection can depend on the accu¬ 
racy of the determination of the continuum shape. In particular, 
we checked if the observed shape can be explained with a mix 
of a stronger reflection edge plus a steeper and lower contin¬ 
uum. However, even if we assume the combination of reflection 
strength and intrinsic spectral shape that minimises the intensity 
of the trough, the normalisation of the Gaussian line is still in¬ 
consistent with the null value at a confidence level of 90%, in 
agreement with the significance of the detection derived from 
the simulations (see below). 

The high energy range where the line is observed, very near 
to the end of the XRT bandpass, could raise concerns of pos¬ 
sible artefacts due to the background. However, the S fN and 
number of source counts collected between 8 and lOkeV (~ 7 
and ~ 56 counts, respectively; last 3 bins in Fig. 2) strongly 
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J ^ ^ ^ .I'o 

Rest Energy (keV) 


Fig. 3. Upper panel: Unfolded XRT data (in black) and average BAT 
spectrum (in blue) fitted including a Gaussian absorption line (E ~ 
7.8 keV, EW-230 eV and cte ~ 190 eV). The continuum (red con¬ 

tinuous line) is modelled with an absorbed plus an unabsorbed power 
law and a reflection component. Lower panel: Relevant residuals, plot¬ 
ted in terms of sigmas with error bars of size one. 

support that the spectral rising that defines the line is real. We 
performed extensive simulations testing the null hypothesis that 
the spectrum is well fitted by a model that does not include the 
7.8 keV absorption feature, as done in Markowitz et al. (2006, 
see also Porquet et al. 2004). Briefly, to take into account the 
observed background and the uncertainty in the continuum, we 
first generated a fake spectrum for a 72 ks exposure assuming 
the best fit found for the continuum. We then fitted this model 
to the fake spectrum, and starting from the new best-fit parame¬ 
ters we re-run a simulation with the same exposure. The baseline 
model has been fitted to this final fake spectrum, and the derived 

has been compared with the minimum value of obtained 
when a narrow (cr = 10 eV) Gaussian component was included. 
We stepped the centroid energy of the absorption line over the 
6.5 - 9keV range in increments of 0.125keV, fitting separately 
each time to derive the lowest value of x^- The whole process 
has been repeated 400 times, and we estimated a 4% probability 
of detecting a similar feature by chance. 

The observed energy of the line suggested by the data is not 
consistent with any of the atomic transitions"^ expected at ener¬ 
gies > 7keV (e.g., Kallman et al. 2004). The most likely ex¬ 
planation is that the centroid of the line is blueshifted, implying 
that the material responsible for the observed feature is outflow¬ 
ing. Prime candidates for the origin of the line are the inner K- 
shell resonances from moderately-highly ionized Fe Ka, as ob¬ 
served in other AGN. Indeed, assuming K/3 absorption by mod¬ 
erately ionized Fe, we would expect to observe also strong Ka 
absorption by the same species (see e.g. the detailed discussion 
in Markowitz et al. 2006). Being conservative, we can identify 
the line with the Fexxvi at £ = 6.966 keV resonant absorption; 
in this case, the observed centroid would indicate a substantial 
blueshifted velocity of (0.11 + 0.03)c. An origin in Fe Ka at a 
lower ionization would imply an even higher blueshift. 

The observed outflow velocity of y out ~ 0.11c would trans¬ 
late in a lower limit on the radial distance, corresponding to 
the escape radius at which the material is able to leave the 
system. When a spherical geometry is assumed, this limit is 

http://physics.nist.gov/PhysRefData/ASD/lines_form.html 


R > Resc - where M is the enclosed mass producing 

the inward gravitational force. Under the reasonable hypothesis 
that this corresponds to the mass of the central BH, assuming 
the estimate reported in the literature (Mbh ~ 1.15 x 10^ Mq, 
see Sect. 2), we have £esc ~ 2.7 x 10'^ cm. This translates to 
a possible launch radius of 10“^ pc, or ~ 80 rs (Schwarzschild 
radii, rs = IGMshIc^) . 

The quality of the available data prevents a more detailed 
spectral analysis, but the ionization level of the Fe responsible 

of the absorption and the observed EW-230 eV (though with 

large uncertainties) would suggest that we are dealing with a 
high-ionization (ionization parameter^ log f ~ 3 ergs cm s“'), 
high-column density (Ah ~ lO^^cm^^) outflow (e.g., Gofford 
etal. 2013). 

4. Summary 

In this paper, we have reported on our 0.3-10 keV observation of 
the Seyfert 1.9 galaxy ISXPS J050819.8+172149, a member of 
the local LIRG pair known as CGCG 468-002. The Swift-XKT 
data have been analysed jointly with the 5w//f-BAT spectrum, 
averaged over 70 months. The continuum is well described by 
an intrinsic power law with photon index F ~ 1.8, absorbed by a 
column of Ah ~ 2.4 x 10^^ cm“^. A reflected component is also 
observed, with a reflection fraction of ~ 1.4, while a weak soft 
scattered component (~ 2% of scattering fraction) can account 
for the observed soft emission. The de-absorbed luminosity is 
T2-i0keV ~ 4 X 10"^^ergs s^‘. 

Our Swift monitoring (the first observation of this source ex¬ 
tending up to energies E ~ 1 - lOkeV), performed in less than 
1 month, provides us a tentative detection (~ 2.1cr significance) 
of an absorption trough at a rest-frame energy of ~ 7.8 keV. 
When the feature is described by a simple Gaussian absorption 
line, its properties (e.g., energy and EW) are consistent with an 
origin in a material moving with a velocity of ~ 0.11c. To our 
knowledge, this would be the first detection with Swift of a pre¬ 
viously unknown high-velocity outflow. 

The low statistics of the data and the high energy of the ob¬ 
served residuals, near the higher boundary of the bandpass of 
XRT, do not allow us to test more physically consistent mod¬ 
els (e.g., grids of photoionized absorbers generated with the XS- 
TAR photoionization code; Kallman et al. 2004). However, if the 
detection is confirmed, the observed EW and the derived veloc¬ 
ity suggest physical parameters typical of an extremely powerful 
outflow, as observed in only a handful of AGN (e.g., the Ul¬ 
tra Luminous Infrared Galaxy/quasar Mrk 231; Eeruglio et al. 
2015). In this case, the kinetic output could match or exceed 
the typical fraction of bolometric luminosities required for AGN 
feedback. 

In fact, 1SXPS J050819.8+172149 could resemble Mrk 231: 
the source studied here is hosted in a star forming merg¬ 
ing system and (possibly) shows evidence of a powerful 
disk wind. However, these characteristics are combined in 
ISXPS J050819.8+172149 with a lower level of activity, both 
of accretion and star formation, than observed in Mrk 231. This 
would make ISXPS J050819.8+172149 quite unique among the 
extremely powerful AGN winds studied so far. Merging systems 
as the one hosting this source are the objects where we expect 
to better observe the interplay between star formation and ac- 

^ The ionisation parameter is defined as ^ = Ljon/nR^, where n is the 
hydrogen number density of the gas (in cm“^) and R is the radial dis¬ 
tance of the absorbing/emitting material from the central source of X- 
ray (in cm), while the ionising luminosity Lion has units of ergs s“*. 
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cretion, since both phenomena can be triggered by galaxy in¬ 
teractions. Indeed, they are the objects were the coexistence of 
disk winds and molecular outflows has been found so far (e.g., 
Mrk 231, Feruglio et al. 2015; IRAS FI 1119H-3257, Tombesi 
et al. 2015). In addition, comparing the accretion and star form¬ 
ing properties reported in Sect. 2, ISXPS J050819.8H-172149 
seems to be one of the few examples of a source that, after a 
recent episode of star formation, is in a transition phase between 
star forming-dominated (HII-LIRG) and accretion-dominated 
(Seyfert-LIRG) state (Alonso-Herrero et al. 2013). The quench¬ 
ing of the star formation can be related (at least partly) to the in¬ 
creasing of the AGN activity, as expected from the co-evolution 
models. An outflow powerful enough to affect the environment 
beyond the SMBH’s gravitational sphere of influence, as the one 
possibly detected in the XRT data, could in principle play a sig¬ 
nificant role in this process. 

Better statistics and higher resolution observations, extend¬ 
ing at energies above ~ lOkeV, are needed in order to confirm 
the presence of the feature, improve the significance of this de¬ 
tection and asses the properties of the associated wind, namely 
the column density and the ionisation state, and then the radial 
location with respect to the central source. The knowledge of 
these parameters would allow us to estimate the mass outflow 
rate and the kinetic power, to be compared with the energetic of 
the accretion. 
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